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We calculate the survival probability of an immobile target surrounded by a sea of uncorrelated 
diffusive or subdiffusive evanescent traps, i.e., traps that disappear in the course of their motion. 
Our calculation is based on a fractional reaction-subdiffusion equation derived from a continuous 
time random walk model of the system. Contrary to an earlier method valid only in one dimension 
(d = 1), the equation is applicable in any Euclidean dimension d and elucidates the interplay between 
anomalous subdiffusive transport, the irreversible evanescence reaction and the dimension in which 
both the traps and the target are embedded. Explicit results for the survival probability of the target 
are obtained for a density p{t) of traps which decays (i) exponentially and (ii) as a power law. In 
the former case, the target has a finite asymptotic survival probability in all integer dimensions, 
whereas in the latter case there are several regimes where the values of the decay exponent for p{t) 
and the anomalous diffusion exponent of the traps determine whether or not the target has a chance 
of eternal survival in one, two and three dimensions. 



I. INTRODUCTION 

Geometric and dynamical constraints imposed by com- 
plex or crowded environments often result in subdiffusive 
behavior, i.e., in sublinear growth of a particle's mean 
squared displacement at long times. However, a com- 
plete description of the underlying transport process at 
a mesoscopic level must go beyond the mean squared 
displacement and involve other properties of experimen- 
tal interest which may be studied via suitable quantifiers 
[U, [1] . This may help one to discriminate between models 
when describing realistic experimental situations where 
subdiffusive (or, more generally, anomalous) transport is 
observed. 

The detailed microscopic subdiffusive transport mech- 
anism is often unknown, and so the literature is popu- 
lated with a number of different models. One popular 
choice to mimic situations of experimental interest is the 
continuous time random walk (CTRW) model 3 with 
a long-tailed waiting time distribution. The CTRW has 
been used successfully as a phenomenological model to 
describe aging effects in systems as diverse as stock mar- 
kets 043 1 charge carrier transport in disordered media 
Q , luminescence quenching in micellar clusters [^[.trans- 
port in porous media [l^, HH , escape problems , and 
morphogen gradient formation |13l4l5l | . 

From a mathematical point of view, the CTRW with a 
long-tailed waiting time distribution and a jump length 
distribution of finite variance is known to be equiva- 
lent to a fractional diffusion equation in the long-time 
limit, that is, a diffusion equation with fractional time 
derivatives rather than ordinary derivatives [l6j . Despite 
the fact that fractional derivatives are non-local integro- 
differential operators, Laplace transform techniques com- 
monly used for the solution of the ordinary diffusion 



equation remain applicable and can be used to tackle 
a wide variety of these problems [17 1 . 

One compelling reason to work with CTRW mod- 
els and the associated fractional equations is that they 
make it possible to include reactive processes. We intro- 
duce this terminology in the broadest sense of including 
particle destruction, creation, binding or transformation 
processes. While the combination of subdiffusion with 
its memory effects and reaction processes is complex, 
at least the CTRW approach offers a way to consider 
them in combination, something that has proved more 
elusive with other approaches. In some fortunate situa- 
tions, the effect of the reactions can be adequately de- 
scribed by suitable boundary conditions imposed upon 
the corresponding fractional diffusion equation (see e.g. 
[ist . Section 4.1 in [13, and references therein); however, 
such situations are rather exceptional, since in general 
the combination of reaction with non-Markovian kinet- 
ics [Tgj leads to non-intuitive fractional equations where 
the parameters describing the chemical kinetics appear 
in a non-universal, model-dependent fashion |2C)| - [24j . In 
particular, heuristic approaches based on fractional equa- 
tions with separate reaction and transport terms such as 
we are accustomed to in ordinary reaction-diffusion prob- 
lems very often lead to unphysical results even in the 
simplest cases of irreversible first-order reactions. 

While fractional reaction-subdiffusion equations have 
been used to investigate a number of different problems 
corresponding to different mesoscopic models and differ- 
ent boundary conditions [lii . [25l - l30i | , many subdiffusive 
versions of classical reaction-diffusion problems 3l| re- 



main unexplored. Thus, one can legitimately claim that 
the field is still in its infancy. One class of problems 
that has attracted considerable interest in recent years 
concerns target search processes driven by (sub)diffusion. 
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Such processes are ubiquitous in nature and include bi- 
nary searches where two objects must meet for a reaction 
or trapping event to occur. In many instances, Smolu- 
chowski's theory of diffusion controlled reactions turns 
out to be a successful tool for the quantitative charac- 
terization of diffusional target search. Examples include 
scavenging reactions [12, Hj] , site location in DNA [sj] , 
ligand binding to sites on macromolecules l35l . predator- 
prey models [36|, luminescence quenching [37[, intermit- 
tent search processes [111, and search processes with re- 
setting to the initial position (39j . to name but a few. 
In this context a key quantity is the so-called survival 
probability of the target, from which the moments of the 
first-passage-time distribution for target annihilation can 
also be straightforwardly computed [13, lllj . 

In recent years the classic diffusional target search 
problem has been generalized to particles that undergo 
anomalous diffusion (42H45| . In this paper we consider a 
related problem, namely, the survival probability of an 
immobile target immersed in a sea of uncorrelated subd- 
iffusive traps that may die "spontaneously" in the course 
of their motion. In other words, there are now two reac- 
tions occurring simultaneously: the disappearance of the 
target and a trap upon encounter with each other, and 
the disappearance of the traps due to some other physical 
process. We term this latter process "spontaneous" as a 
way to recall that it is not induced by collision with the 
target. The spontaneous evanescence process may for in- 
stance be triggered by particle scavengers in the system, 
but for practical purposes any process that turns off the 
interaction between a trap and the target can also be 
thought of as an evanescence or death process. 

A solution to this problem in dimension d = 1 was 
given in Rcf. [i^ using a functional method first devel- 
oped by Bray et al. [43] for the diffusive case. Here we 
approach the problem from a different point of view that 
allows us to also obtain results in higher dimensions. In 
particular, we make use of a recently derived reaction- 
subdiffusion equation obtained from a mesoscopic CTRW 
model with a long-tailed waiting time distribution and a 
superimposed reactive process. The elimination of the 
fixed target is incorporated as a boundary condition, 
while the decay mechanism of the traps as they move 
subdiffusively is modeled by a phenomenological choice 
of a monotonically decaying functional form for the trap 
density p{t). 

The paper is organized as follows. In Scc.|ll]we briefly 
recall the results for the survival probability of a target 
in a sea of non-evancscent traps. Our new general results 
for the survival probability of the target when the traps 
are evanescent are presented in Sec. lIIIl In Sec. lIVI we im- 
plement these results for particular forms of evanescence, 
namely, exponential and power law. We conclude with a 
short summary in Sec. fVl 



II. RECAP OF RESULTS FOR 
NON-EVANESCENT TRAPS 

We consider a statistical ensemble of systems each of 
which is composed of a fixed hyperspherical target of ra- 
dius R located at the center r = of a large d-dimensional 
volume V. Each volume V initially contains A^o randomly 
distributed non-interacting point traps. At time t = the 
traps begin to move subdiffusively. If any of them hits 
the surface of the target, both the target and the trap are 
instantaneously annihilated (fully absorbing case). Our 
goal is to compute the survival probability of the target 
at time t, i.e., the probability that no traps have col- 
lided with the target up to this time. In this section, we 
briefly recall previous results obtained in Ref. [42| when 
the traps are not subject to spontaneous evanescence. In 
the next section we use these results to obtain the solu- 
tion when the subdiffusive traps evanesce. 

In the absence of evanescence, the motion of each trap 
is dictated by the fractional diffusion equation 

= V? -(r, t|ro; 0), < 7 < 1, 

(1) 

where w{r,t\ro;0) is the probability density of finding 
the trap at location r at time t if it started at position 
To at t = 0, K-y is the anomalous diffusion coefficient, 
and Vr stands for the Laplacian operator with respect 
to the position r. The operator oT>l~'^ is defined via the 
equation 

A^n {oVl-\fit)} = u'-'Ct^u {.fit)} , (2) 

where £t^„{/(i)} = f{u) = e~^*- f{t) dt denotes the 
Laplace transform (the function and its Laplace trans- 
form are clearly distinguished by the argument and so 
we use the same designation for both). Instead of o1^t~'^ 
we employ the more commonly used Riemann-Liouville 
fractional derivative qDI^^ defined as follows [48| : 

oD]-'' f{r, t) = f dt' /^':!? ■ (3) 

* ' r(7) dt Jo {t - t'y-'i ^ ' 

Strictly speaking, the two operators are not identical. 
However, they are equivalent for functions which are suf- 
ficiently regular at i = 0, as are all the functions we en- 
counter in this problem. The propagator solution of the 
fractional diffusion equation ([1]) yields a mean squared 
displacement with the long time behavior (r^) ~ AT-yP re- 
sulting in subdiffusive behavior when 7 is less than unity. 

Let QT{t;R) denote the ensemble averaged survival 
probability of the target in a sea of randomly distributed 
uncorrelated traps. This quantity can be obtained from 
the survival probability of the target in the presence of a 
single trap starting at location rg, Qi^Tij^o,t] R). We 
shall focus on the thermodynamic limit, i.e., we take 
A^o cxD, — > 00 while keeping a fixed global initial 
trap density po — VmYMa,v-^oo Nq/V . In this limit one 
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QT{t]R) 



hni 

No.V^oo 



j7 I Qi.T{ro,t;R)dro 

^ Jra>R 



No 



(4) 



where the integration is carried out over the volume that 
is exterior to the target. We have introduced the auxil- 
iary quantity 

<j{t; R)^4d [ [1 - Qi,T(ro, i; R)] dro 

^ Jro>B, 



p^Rd- 



\nQT{t-R). 



(5) 



Note that the survival probability (5i_r(ro, t; i?) of 
the target is identical with the survival probability 
Qi{Y(),t\R) of the trap, as we have assumed that both 
the target and the trap disappear instantaneously upon 
encounter, i.e., there is only one decay channel for both 
particles. Note also that this is no longer the case when 
the traps undergo spontaneous evanescence. In that case 
one has Qi_t > Qi (sec next section). 

In order to compute Qi{YQ^t] R), we must first define 
the relevant boundary value problem by complementing 
Eq. ([1]) with the deterministic initial condition. 



-u;(r,0|ro;0) = 6{v - tq). 



and the boundary conditions 



w{R,t\vo-Q) = 0, 
lim w(r, t|ro; 0) — 0. 



(6) 



(7) 
(8) 



The boundary condition ([7]) reflects the fully absorbing 
nature of the target, which prevents the trap from being 
found on the target surface or inside the target. The solu- 
tion ^(r, t|ro; 0) is related to (5i(ro, t; R) via the integral 
relation 



Qi(ro,i;-R)= w{r,t\ro;0) dr. 



(9) 



The spherical symmetry of the target means that 
Qi(rQ,t; R) only depends on the initial distance tq of 
the trap to the target. For this reason, we shall drop the 
subindex of ro and from here on use the simpler notation 
Qi{r,t; R). Taking into account Eq. ©, the boundary 
value problem stated directly in terms of Qi{r, t\ R) then 



dQi{r,t;R) 



at 

Qi(r,0;i?) = l, 
Qi{R,t;R) = 0, 
lim Qi{r,t]R) = 1. 



(10a) 

(10b) 
(10c) 
(lOd) 



The second equation in this set corresponds to the initial 
condition and is self-explanatory, while the third one is 



a boundary condition which reflects yet again the fully 
absorbing nature of the target. The last equation states 
that a trap which is "pushed" infinitely far away from 
the target will survive forever, as its only decay channel 
is provided by the interaction with the target. 

The above problem can be solved exactly in Laplace 
space (i^ . 



-I Kd/2-i 

Ka/2-1 i^/RhrfjK. 



uQi{r,u;R) ^ - 



(11) 

where i^d/2-i(") is a modified Bessel function of the 
second kind. Alternatively, the solution Qi{r,u; R\j) 
for 7 < 1 can be found from the corresponding so- 
lution for normal diffusion (7 = 1) by means of the 
"time-expanding transformation" [l^H^I associated with 
the so-called subordination principle: uQi{r,u; R\j) = 
u'>Qiir, uT;i?|7 = l). 

For d = 1 and d = 3 the Bessel functions can be ex- 
pressed in terms of exponentials, and explicit exact solu- 
tions are available for arbitrary times t. In other dimen- 
sions simple expressions are only available at long times. 
Setting Qi,T = Qi in Eq. ([3]), Tauberian theorems can 
be used to find the long-time behavior in the multiple 
trap problem [i^ : 



a{t;R) oc 




d=l 
d= 2, 
d>3 



(12) 



2U/7 



where = {AK^/R^) 

Thus, the survival probability of the target goes to zero 
in all dimensions d. This result is in strong contrast with 
the single-trap problem, since in that case the probability 
that the random walk of the trap never intersects the 
target becomes non-zero as soon as d > 3. 



III. SURVIVAL PROBABILITY FOR 
EVANESCENT TRAPS: GENERAL 
EXPRESSIONS 

The behavior of the survival probability of the target 
changes completely if the traps disappear in the course of 
their motion. We assume a spontaneous evanescence pro- 
cess, specifically, that the decrease of the global trap den- 
sity p{t) is described by the following differential equa- 
tion: 



-A(t)p(t), 



(13) 



where \{t) > is a rate coefficient which is in gen- 
eral time dependent. The solution po exp ^— Jp* \{t') dt'^ 

yields a decaying density of surviving traps which de- 
scribes the time evolution of the trap density in the ab- 
sence of the fully absorbing target. The case A = constant 
leads to an exponentially decaying density. 
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Our main goal is to compute the survival probability 
Q^{t; R) of the target at time t (we use survival prob- 
abilities with a star to distinguish these quantities from 
their counterparts in the absence of trap evanescence). 
We follow the strategy of the previous section, namely, 
to derive the solution from the single-trap case. 

We wish to combine the effects of Eq. and Eq. ([1]) . 
One might be tempted to proceed as in the case of or- 
dinary diffusion and simply construct some superposi- 
tion of transport and reaction terms. However, a careful 
analysis shows that this is incorrect. Instead, a rigorous 

I 



derivation starting at the level of the CTRW shows that 
the correct equation is [l^ 



dw{r,t\rQ,0) _ p{t) 



dt 



Po 



1-7 Po V72 



Pit) 

PW 



V>(r,t|ro,0) 



(14) 



It is straightforward to show that the survival probability 
Ql{ro,t;R) = Q\{rQ,t;R) = / ?ii(r, t|ro; 0) dr of the trap 
then obeys the equation 



dQl{r,t-R) 
dt 



Po 



-7 ( _P^^2 



)l{r,t;R) + ^^Ql{r,t;R), 



(15a) 



where we have again dropped the subindcx from ro for 
notational convenience. The above equation must now 
be complemented with the conditions 



Qlir,0;R) = l, 
Ql{R,t;R) = 0, 

Pit) 



lim Ql{r,t;R) 



PO 



(15b) 
(15c) 

(15d) 



Note that the last equation is substantially different from 
the corresponding one in the absence of the evanescence 
process. Indeed, even the survival probability of a trap 
which is at an infinite distance from the target decays in 
time because of the evanescence reaction, and this prob- 
ability is equal to the ratio of the global trap density at 
time t in the absence of the target and the initial density 
Po- 

At this stage, one can easily check that if one per- 
forms the transformation QUr, t; R) = [pit)/ po]Q{i', t; R) 
in Eqs. ([15]), the resulting set of equations for the trans- 
formed function Q is identical with the boundary value 
problem pO)) for Qi. Because of the uniqueness of the 
solution, we thus conclude Q = Qi, and further. 



QUr,t;R) = ^Qi{r,t;R). 
Po 



(16) 



This expression is intuitively clear: it simply states that 
the probability that up to time t the trap has neither 
evanesced (one decay channel) nor hit the target (another 
decay channel) is equal to the product of the probability 
pit) / Po that the trap has not evanesced, and the condi- 
tional probability that it has not hit the target given that 
it has not previously evanesced. The latter probability 
is precisely the survival probability of the trap when no 
evanescence is at play. 

The next step in our route to the solution for the multi- 
trap problem is to derive a relation between ^(r, t; R) 
and Ql{r,t; R). As already anticipated in the previ- 
ous section, the evanescence of the trap implies that the 



I 

survival probability of the target and the trap are no 
longer the same. Let t' be a time in the interval (0,t). 
In the single-trap problem, the infinitesimal probabil- 
ity {d[l - Qlrp{r,t';R)]/dt'}dt' that the target is an- 
nihilated by collision with the trap during the interval 
{t', t' + dt') is the product of two factors: i) the probabil- 
ity {d [1 - Qiir, t'; R)]/dt'}dt' that the trap coUides with 
the target during the time interval (t' , t' + dt') given that 
it has not previously evanesced, and ii) the probability 
p{t')/po that up to time t' the trajectory of that trap is 
not interrupted by an evanescence event. Thus, 



d , 



Ql^{r,t';R)]dt' 

Po dt' 



i{r,t';R)]dt'. 

(17) 

We next implement a number of steps [integrate this 
equation, integrate by parts, and use Eq. (|16p ] to ob- 
tain 



Qyir,t;R) = QUr,t;R) 



' Ql{r,t';R)^ dt', 
pit') 



(18) 

which quantifies the difference between Q* rpir, t\ R) and 
Q\ir,t-R). 

Having clarified the relation between Ql{r,t;R) 
and Qi{r,t]R) (survival probability of a single non- 
evanescent trap) and between Qlj,{r,t;R) (survival 
probability of the target in the presence of a single 
evanescent trap) and Q\{r,t;R) (survival probability of 
a single evanescent trap) , we are now ready to tackle the 
multiple trap problem by proceeding as in the previous 
section, i.e., by using the statistical independence of the 
traps. Equations (|4|) (for the survival probability of the 
target in the presence of a collection of traps) and ^ are 
now respectively replaced with 



Qxit] R) 



lim 



Qlj.{r,t;R) dr 



r>B. 



No 



:exp{-poi?V(t,i?)}, 



(19) 
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and 



= [l-QlAr,t;R)]dr 

^ Jr>R 

^ ■\nQ*T{t;R). 



(20) 
(21) 



Taking the derivative of Eq. ([20|) with respect to time 
and using Eq. ([T8| we get 



da*{t;R) 

di 



1 

R^ 



r>R 



dt p{t) 



dr. 



(22) 



Next we use Eq. (|15ap and the relation ([TS]) in the right 
hand side of Eq. ([^^ and apply Gauss' theorem to change 
the volume integral to a surface integral. This allows us 
to write 



da*{t-R) 
di 



S, 



Kj P{t) dQi{r,t;R) 
R po 



dr 



(23) 



r=R 



where Sd = 27r'*/^/r(c?/2) denotes the surface of a d- 
dimcnsional hypcrsphcre of unit radius. Finally, inte- 



grating from to t and using the condition cr*(0; i?) = 
[Q* t(^' ^) ~ 1] obtain the general formula 



a*{t;R) = Sd^ 

^ JO 

Po 



oD 



dQi{r,t';R) 



dr 



r=R. 



(24) 



which extends the result for non-evanescent traps ob- 
tained in [iH to the case of evanescent traps. More 
specifically, we have again related the logarithm of the 
survival probability of the target, now in the presence of 
a collection of evanescent traps, to the survival probabil- 
ity of the target in the presence of a single non-evanescent 
trap. If the trap density p{t) decays sufficiently rapidly, 
then Qlxi^ ~^ oo,R) > and a*{t — >■ oo;i?) < co, and 
it is possible to conclude from Eq. that 



a*{oo;R) -a*(t;R) 
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dQiir,t';R) 



dr 



r=R. 



PO 



dt'. 



(25) 



To continue from here we need to specify the dimension- 
ality explicitly, and also the explicit form of the trap 
density as a function of time. First we discuss the di- 
mensionality. The case d = 1 has been dealt with in 
[46j ; suffice it to say that we recover the result obtained 
therein: 



r 



With this result we get from Eq. 



(T*{t-R) 



7 



Por(7/2) 



1 2 

(y*{t-R) = — 



poRTi'j) Jo 



p{t')t'^/^-^dt' 
p{t')t'^-^ dt' . (30) 



poR r(7/2) 



/ p{i') t'^^"^ ^ dt' , (26) Here we have used the fact that the fractional derivative 

JQ r-i :„ u..-<- : 



where we have made use of the explicit one-dimensional 
form of Qi in terms of the Fox iJ- function 



Qi(r,i;i?) = l-i/io 



r-R 


(1,7/2)" 




(0,1) 



(27) 



for the survival probability of a non-evanescent trap. For 
d = 3, and exact expression for the survival probability 
valid for arbitrary t can also be obtained. This is a new 
result so we proceed in more detail. For the single-trap 
problem without evanescence we have 



of C, a constant, is not zero but is instead qD] C 

cp-Vr(7). 

In contrast with the d = \ and the d = 3 cases, no 
explicit solution in a simple integral form similar to that 
of Eqs. (gH]) and is available for a*{t;R) when d 
2. However, we can write explicit expressions for the 
approach of the survival probability to its final value. 
From the result for Qi given in [42|, one finds 



dQi{r,t;R) 



dr 



1 



r=R 



R 7 In(a-yi) 



(31) 



Qi(r,t;i?) = 1 

leading to 

dQi{r,t;R) 



dr 



r=R 



1 

R 



r- R 


(1,7/2)' 




(0,1) 


1 


1 



r(i-7/2) 



(28) 



(29) 



The fractional derivative of this expression is best com- 
puted in Laplace space [l^j. Transforming the resulting 
expression back into the time domain we find that 
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[l/ln(a^t)]^t^-V[r(7)ln(a^i)]- 



(32) 



Let us further assume that the trap density p(t) de- 
cays sufficiently rapidly to ensure that a*{oo; R) is finite. 
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From Eq. ([25]) wc find that for d = 2, 



a* {oo; R) ~ a* {t;R) 



poR^T{-f + 1) Jt ln{a^t') 



We next implement our general results for particular 
forms of the decay of the trap density. 



IV. IMPLEMENTATION FOR PARTICULAR 
TRAP DECAY FUNCTIONS 

A. Exponentially decaying trap density 

We first consider an exponentially decaying trap den- 
sity, p{t) = poexp(— At), where the characteristic time 
scale of the decay is given by r = A^^. This would rep- 
resent a unimolecular decay if this were the only decay 
channel, which in the presence of the target it is not. 

I 



However, for a single particle, or for the first of many, 
this is still the scenario because we are looking at the 
decay only up to the time that this second channel first 
affects the trap density. The focus of our interest is in the 
final value of the survival probability of the target and 
the long time approach to it, but for completeness we also 
give a general expression for the early time behavior. 

For d = 1 and d = 3 wc can directly insert the expo- 
nential decay of the density in Eqs. ((26|) and (pO)) . and 
perform the integrals. Alternatively, wc note that for the 
exponentially decaying trap density Eq. leads to the 
relation 

Qt,r(r,t;i?)-e-^*Qi(r,t;i?) 

+A / e-^*'Qi{r,t';R)dt'. (34) 

Next, from Eq. (fTT]) and the Laplace transform of 
Eq. ([M]) it then follows that 



uQl rp{r,u;R) = 1 



f- 

\R 



y 



Ka,2-iWR^{u + XVlK^)' 



(35) 



r 



The Laplace transform of Eq. 
leads to 



(PHI then immediately 



(T*(w; R) 



K^,2WR^{u + \)'</K^) 



uR{u + A)^/2 Kd/2-iWR''{u + \VlK^) ' 



For odd- valued d (but not for even- valued d) the modified 
Bessel functions of the second kind can be expressed more 
simply in standard power series expansions. For d = 1 we 
can follow either path (direct integration or simplification 
and inversion of the modified Bessel function) to obtain 



IP 

a*{t;R)^^ 



1 



r(7/2,At) 

r(7/2) 



(37) 



where £^ 



is a characteristic length scale 



associated with the distance covered by a trap during its 
mean survival time when no tar get is present. Eq. ([57)) 
is equivalent to Eq. (15) in Rcf. [4Q|. For d = 3 we can 
again follow either route to the solution and find 



a*{t;R) = 47r 



R 



1 - 



£2 



r(7/2,At) 

r(7/2) 
r(7,At) 
r(7) 



1 - 



resulting in a smaller survival probability than in d = 
1. As noted already, there are no simple, closed- form 
solutions valid for arbitrary times for d = 2, although we 
are able to extract some limiting behaviors for this case 
as well (see below). 



The computation of the final value of a* {t; R) in arbi- 
trary integer dimension is readily obtained from ([36]) by 
means of the final value theorem for the Laplace trans- 
form: 



(36) fT*(oo;i?) 



lim u (t" 

tt->0 



iu;R) 



f 



Kd/2{R/t^ 



RKa/2-,{R/e^) 



(39) 



Before further evaluation, we note as an aside that the 
above non-zero survival probability implies an infinite 
mean survival time of the target in any dimension d. In 
contrast, if the traps do not evanesce, the mean lifetime 
of the target is finite [l^]. For odd dimensions, Eq. ([M)) 
can be conveniently rewritten as rational functions of the 
argument R/i^. For example : (T*(oo;i?) = 2i^/R for 
d ^ 1 and cr* (oo; R) = (M^/R) (1 + l^/R) for d = 3. 

Next we explicitly present the results for the ap- 
proach to the final value by exhibiting the difference 
cr*(oo; R) - a*{t; R) at long times. For d = 1 and d = 3 
this respectively follows from Eqs. ([57)) and ([55)) . whereas 
from d = 2 the long time behavior can be inferred from 
Eq. We find: 



1 



(38) a* {oo; R)—(T* {t; R) oc < 



In ^{a-yt)^^ 



d = 3. 

(40) 

This asymptotic behavior also holds in the case of normal 
diffusion (7 = 1). One can see that the decay of the sur- 
vival probability to the final state prescribed by Eq. ([^0)) 
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becomes faster as one goes from one to two dimensions 
and from two to three dimensions. It is also straightfor- 
ward to show that the long time behavior prescribed by 
PO)) for d = 3 remains vaUd for d > 3 (the prefactor, 
however, depends on d). 

The hmit A — > (no evanescence) turns out to be 
singular. Indeed, in the absence of evanescence, a*{t] R) 
tends to infinity as given by Eq. which is different 

from the result obtained when taking the limit A — ^ 
in Eq. pO)) . We conclude that the evanescence reaction 
completely changes the physics of the problem, affecting 
both the steady state and the decay form of the survival 
probability. 

Finally, the short time behavior {t X^^) is straight- 
forward to obtain via a Tauberian theorem applied to the 
large u (w 3> A) limit of Eq. (fMl) : 



a*{u;R) 



a*{t;R) 



i?r(l + 7/2) 



(41) 

As one might have guessed, the short time result is inde- 
pendent of A, i.e., the effect of the evanescence reaction 
is still negligible in this regime. 



B. Power law decay of the trap density 

We next turn to the case of a power-law decaying den- 
sity, that is. 



The behavior of the integral on the right hand side de- 
pends on the relative values of /3 and 7. We consider 
three different cases: 

Case 1 (/3 > 7). From the asymptotic long-time be- 
havior pip and the expression for the fractional deriva- 
tive of the inverse logarithm ([5^ one finds 



dQi{r,t'-R) 



dr 



r=R 



i?r(7 1) In(a^i') ' 



(45) 



Under the assumption that t is large enough, we now 
split the interval of integration [0, t] into two subintervals, 
Ii = [0,ic] and I2 = [tcji)i where tc is chosen sufficiently 
large so as to ensure that the approximation (P5|) holds 
over the full extent of l2. Hence one has 



cr*(i;i?) 



r(7-|-l)i?2 J^^ ln{a^t' 



dt' 



(46) 



where we have used the long-time approximation 
p{t')/po ~ {i' and C represents the integral from 
to In this case one can easily check via partial inte- 
gration that the integral on the right hand side of Eq. p6)) 
remains finite as i — > 00. Hence cr*(oo; R) is finite and the 
target has a non-zero chance of eternal survival. Using 
the explicit form of p{t') in Eq. we find 



a*{oo;R) - a*{t;R) 



(/3-7)r(7 + l)i?2 \x,{a^ty 



(47) 



Pit) 



Po 



{l + t/T)P-' 



/? > 0. 



(42) 



This choice corresponds to a time dependent rate con- 
stant A(t), which can be used to capture the essential 
features of complex higher-order kinetics by means of the 
linear differential equation (|13p and a proper choice of the 
exponent f3. Interestingly, the survival probability of the 
target in this case depends not only on dimensionality 
but on the relative values of the power law decay expo- 
nent /3 and the subdiffusion exponent 7. 

The behavior when d = 1 follows directly from Eq. 
I, and can be summarized as follows [46| : 



. r(/3-7/2) 

■Ht/r), 



a*{t;R) 



i?,r(7/2) 



/3>7/2 
/3 = 7/2 
{t/rr/^-P, /3 < 7/2 

(43) 

where is defined in a way similar to the exponential 
case, i.e. , -^-y 



fl(7/2-/3)r(7/2) 



(X^r^)i/2. 



Next we consider the two-dimensional system. Our 
starting equation is with d = 2: 



a*{t;R)^2TT^-^ 

n 70 



x{l-t'/T)-^ dt' 



dQi{r,t';R) 



dr 



r=R. 



(44) 



Unfortunately, it does not seem possible to find an ex- 
plicit exact expression for ct*(oo; R) due to the lack of an 
exact expression for (5i(r, t'; R) valid for the whole time 
domain. 

Case 2 (/3 = 7). In this marginal case the target also 
disappears eventually, but the approach to the empty 
state has a different analytic dependence, as Eq. p6)) 
now leads to 



r(7 + l)i?2 



ln(ln(a7<)). 



(48) 



Thus, the target is eventually killed with certainty, in 
agreement with the result given in Ref. for the special 
case /3 = 7 = 1 (normal diffusive traps). 

Case 3 (/3 < 7). Clearly, in this case the constant C of 
Eq. (Uni) becomes negligible at sufficiently long times and 
the behavior of <J*{t\ R) is dominated by the integral on 
the right hand side. Using partial integration one easily 
sees that the leading long-time behavior is given by 



a*{t;R) 



(7-/3)r(7-Hl)i?2 In(a^t)' 



(49) 



Hence when /3 < 7 the target eventually disappears with 
certainty. 

We next discuss the three-dimensional case. In this 
case, the integrals in Eq. (PH)) can be carried out exactly. 
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and one finds 

a*{t;R) 



J^^t/(r+t)(7/2,/3-7/2) 



(50) 



where 



is the incomplete Beta function [SlJ. The long-time be- 
havior of a*{t;R) again depends on the relative values 
of P and 7. The analysis is carried out along lines sim- 
ilar to those presented in Ref. [i^ for d = 1. We shall 
distinguish three different cases. 



/•a 

Ba{z,w) = dt with Rc(z) > 0. (51) Case 1 (/3 > 7). We can rewrite (j*{t;R) as 

Jo 



<^*it; R) = 3(7/2, /3 - 7/2)/,/(,+,)(7/2, /3 - 7/2) + J^^B{^^ P - l)It/(r+t){l. P - l)- 



(52) 



Here B{z, w) is the Beta function (where the requirement 
Re(z) > and Re(w;) > places us in the "Case 1" 
regime), and Ix{z,w) is the regularized incomplete Beta 
function as defined in Sec. 6.6.2 (pg. 263) of Ref. fsif . 
Using the property 6.6.3 in [5l| we can set Ixia, b) = 
1 — Ii-x{b,a). Applying the relation 26.5.5 in [5l|, and 
making use of the relation between the Beta function and 
the Gamma function, we arrive at the asymptotic result 



Thus, in two and three dimensions the target has a 
finite probability of surviving forever only for /3 > 7. 
For comparison, in the one-dimensional case it was found 
that the target has a chance of eternal survival only when 
/3 > 7/2 [cf. Eq. (HSl)]. We thus see that the interplay 
between subdiffusivc transport and the evanescence re- 
action determines, also in dimensions higher than one, 
whether the target can ultimately survive. 



a*{t;R) - (T*{oo;R) - 



(t/r) 



R (/3-7)r(7) 



(53) 



SUMMARY AND OUTLOOK 



with 



a*{po\R) = ^-kL 



r(/3-7/2) , 7r^2r(/^_^) 



r(/3) 



R r(/3) 



(54) 



leading to a non-zero survival probability Q^{oo; R) — 
exp [-po-R^f *(oo; R)]. 

Case 2 (/3 = 7). In this case the incomplete Beta 
function in the term proportional to in Eq. (|50p can be 
rewritten as a hypergeometric function and consequently 
for long times this term can be approximated by 



i?r(7 + i) \t 



2^1(7,7,7+1, -t/r) 



i?r(7) 



(55) 

On the other hand, the term proportional to £j goes 
to a constant for long times, as can be seen using the 
same expansion as the one used for the /3 > 7 case. 
Hence, the survival probability Q^{t; R) vanishes as 
{t/TyP"^^V^^i)^ that is, a*{t;R) oc ln(t/r). 

Case 3 (/? < 7). In this case the term propor- 
tional to can easily be seen to behave as [7r£^/i?(7 — 
f3)T{'-f)]{t/T)''~^ by performing a straightforward asymp- 
totic analysis of the corresponding integral. On the other 
hand the £.y term is negligible compared to the £'i^ term. 
This results in a stretched exponential decay to zero, 
i.e., Q^{t;R) oc exp (-Ci''-'^) with C > 0, that is, 
a*{t;R) oc T"^. 



We have presented a particular application of a re- 
cently derived fractional reaction-subdiffusion equation, 
namely, the study of the behavior of the survival proba- 
bility of an immobile target surrounded by a sea of non- 
interacting diffusive or subdiffusive point traps subject to 
an evanescence reaction. The evanescence reaction is as- 
sumed to take place independently of the CTRW jumps 
performed by the traps, as opposed to a recently intro- 
duced model where disappearance takes place at the time 
of each jump js^ . 

The problem considered in this paper is only one of a 
family of many possible boundary value problems which 
may be dealt with using our equations. However, this 
particular choice may be of interest in a number of ex- 
perimental situations, e.g., radical recombination kinetics 
in the presence of added scavenger molecules [s^ also re- 
sponsible for the disappearance of radicals. As far as we 
know, the interplay between the scavenging reaction and 
possible memory effects arising in some environments re- 
mains unexplored. 

We focused on the case of exponential evanescence and 
power law evanescence, extending previous results appli- 
cable only to the one-dimensional case. In particular, our 
results also hold for the normal diffusion case (7 = 1). 
The presence of the evanescence reaction was found to 
completely modify the physics of the problem, both at 
the level of the steady state and the decay of the sur- 
vival probability to a finite steady state or to zero. More 
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specifically, with an exponentially decaying trap density 
p{t) = pqC"'*'* (with A > 0), we find that there is a fi- 
nite survival probability of the target in all dimensions 
because the traps die sufficiently quickly in their search 
of the target. By way of contrast, when the traps do 
not evanesce the target has a zero survival probability in 
all dimensions. The long-time approach toward the final 
value of the survival probability turns out to be more 
complex than in the case of non-evanescent traps, and 
in the subdiffusive case 7 < 1 it involves powers of t as 
well as exponential factors e""^* (with a logarithmic cor- 
rection in d = 2). On the other hand, when the density 
decays as a power law, p{t) oc with /3 > 0, the behav- 
ior depends on the relative values of /3 and the anoma- 
lous diffusion exponent 7 of the traps. In one dimension, 
the target has a finite asymptotic survival probability if 
/3 > 7/2, whereas in two and three dimensions the target 
only has a finite chance of eternal survival when /3 > 7. 

A natural extension of this work would allow normal 
diffusive or subdiffusive target motion (the case of nor- 
mal diffusive target and normal diffusive evanescent traps 
has been considered in Ref. [1^). Note, however, that in 
such a case the respective distances between the target 



and the traps would no longer evolve as independent vari- 
ables, implying that our asymptotically exact approach 
would not work in its present form. Nonetheless, approx- 
imations based on the fact that at long times the domi- 
nant contribution to the survival probability comes from 
the subset of trajectories where the target remains im- 
mobile [13, [13 could prove useful to tackle the problem. 
Ultimately, this behavior finds its roots in what has been 
termed the "Pascal principle" in the literature [56l - [58| . 
i.e., a target placed in a symmetric initial distribution 
of traps survives longer on average if it stays still rather 
than if it moves. 
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